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Abstract 

We study the two-dimensional renormalization-group flow induced by perturbations 
that reduce the global symmetry of the 0(3) cr-model to the discrete symmetries of 
Platonic solids. We estimate the value of the correlation length at which differences in 
the behaviour of the various models should be expected. For the icosahedron model 
with nearest-neighbor interactions, we find £ > 200. We provide an explanation for the 
recent numerical results of Patrascioiu and Seiler and of Hasenfratz and Niedermayer. 
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1 Introduction 



Quantum field theories with a non-Abelian continuous symmetry group of invariance 
play a major role both in particle and in condensed-matter physics. Two cases are of 
particular interest: two-dimensional spin models with non-Abelian global symmetry 
group, and four-dimensional gauge theories with local non-Abelian gauge invariance. 
According to the common wisdom these two cases share the peculiar feature of asymp- 
totic freedom (AF). 

In the lattice formulation it is quite easy to replace the continuous group by one 
of its discrete subgroups. In this case, due to the discreteness of the group, the action 
has a finite gap and at least a freezing transition is expected.^] Nonetheless, at large 
enough temperature one expects only small differences and indeed, finite subgroups 
have been used in Monte Carlo updates to simulate continuous groups ||, 

When the symmetry group is Abelian it may happen that the discrete symmetry is 
enlarged to a continuous one. The massless intermediate phase of Z{N) models with 
N > 4 is the same of the 0(2) model for low enough temperature [||, ||]. 

It has also been suggested that a similar phenomenon occurs in non-Abelian models. 
Patrascioiu and Seiler [||, [?], ||, [9| have often criticized the conventional wisdom on AF for 
the continuous group and proposed an alternative scenario in which a low-temperature 
massless phase appears. If this possibility really happens it is plausible to accept the 
idea that, for example, in d = 2 the 0(3) model is in the same universality class of the 
dodecahedron spin model. 

In order to test this conjecture several large-scale simulations have been performed 



|1C, 11, gj. In particular, the finite-size scaling curve for the second- moment correlation 
length measured in ]i~2| , |i~3| and the renormalized coupling for the 0(3) model have been 
compared with the results obtained for discrete spin models with different discrete 
subgroups and nearest-neighbor interactions. It was found that the icosahedron model 
and the 0(3) model are practically indistinguishable at present-day values of £. This 
is not totally surprising: after all, for (3 small enough one expects only tiny differences 
since in the presence of large fluctuations the discreteness of the spin space should not 
play an important role (this was indeed the motivation of Refs. |Q, [|). However, what 
is more surprising is that the discrepancy seems to decrease as (3 increases (see the 
results of Ref. for the renormalized zero- momentum four-point coupling), while 
naively one would have expected the opposite. 

These numerical results have been interpreted as evidence that the 0(3) and the 
icosahedron model have the same continuum limit. In this paper, we will show that 
this conclusion is in contrast with the common theoretical understanding of the 0(3) 
model: If the continuum limit of the 0(3) a model is correctly described by the pertur- 
bative renormalization group (RG), then the icosahedron and the 0(3) model belong 
to different universality classes, contrary to what suggested in Refs. [|To| , pjj[ . 

For this purpose, we study the effect of perturbations that break the 0(3) symmetry 
down to a discrete subgroup and show that any such perturbation is a relevant pertur- 
bation that modifies the universal behavior. A similar analysis for the cubic symmetry 
was performed by Pelcovits and Nelson p3] and in the context of the XY model in Ref. 



. We then try to provide an explanation to the numerical data. We analyze a model 



1 In some models, for instance in some Z(N) models, between the ordered and the disordered phase, there 
may also be an intermediate massless phase [jjj . 
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with icosahedral symmetry that interpolates between the nearest-neighbor 0(3) model 
(he = 0) and the standard icosahedron model (h§ = +00), using perturbation theory. 
When /i6 is small, we find for ^ < 200 a behavior that is similar to that observed 
numerically: the difference between the icosahedral and the 0(3) model decreases as j3 
increases. However, this apparent convergence is misleading, because for £ > 200 the 
difference between the two models increases as (3 increases, with the result that the two 
models do not lie in the same universality class. On the basis of these results we predict 
that the difference between the standard icosahedron model considered in jnj [ll], |8| 
and the 0(3) model should probably become visible only for £ > 200. 

We stress that all our arguments assume the validity of the perturbative RG for 
the 0(3) model. If this is not correct, it is possible that the icosahedron and the 0(3) 
models have the same universal behavior as conjectured by Patrascioiu and Seiler. 



2 Renormalization-group analysis of perturba- 
tions with discrete symmetry 

In this Section we want to perform a RG analysis of the discrete models. For this 
purpose we consider the lattice Hamiltonian 



latt 



(3H 



where cr T is a unit vector in H 



(1) 



' X+fl 



(2) 



and [i = 1, 2 are the positive directions on a square lattice with lattice spacing a = 1. 
Here, I n (a x ) is a polynomial of cr x that is invariant under the action of the discrete 
group and belongs to an irreducible representation of the 0(3) group, i.e. it has a 
well-defined 0(3) spin n. Such polynomials can be found for all discrete subgroups 
that correspond to Platonic solids. The lowest-rank tensors for each subgroup are 



obtained in App. A. 2. For the tetrahedron, the cube (or the octahedron), and for the 
icosahedron (or the dodecahedron), the lowest-rank invariant tensors have spin n = 3, 
4, and 6. Explicitly, we have for the three cases respectively 



/ 3 
h 




-i3xl=-(Y 3 , 2 - y 3 ,-2), 
35 




Y^ + J—{Y AA + Y A> ^) 

y 6 , q - ^(y 6 ,5-y 6 ,-5: 



Here, we have introduced the spherical harmonics 



(3) 



Y ljT , 



'21 + 1(1 -m) 
47r (I + m) 



P;, m (cos6>) e f 



(4) 
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where Pi )Tn (x) are the associated Legendre functions, and we have written the three- 
dimensional spin <r in terms of polar coordinates 9 and <fi. These tensors are not unique, 
since the normalization is arbitrary and they are defined modulo a rotation. f\ For the 
purpose of showing that discrete-subgroup perturbations are relevant, the value of n 
plays no role and one can use any invariant tensor. One could also consider tensors that 
do not have definite spin and are sums of invariant 0(3)-irreducible tensors. This would 
only make the discussion more cumbersome without changing the physical results. 
The partition function is defined as usual by 



(5) 



The Hamiltonian (|l]) interpolates between the 0(3) model (h n = 0) and the discrete- 
symmetry model (\h n \ = oo).H 

If, as claimed in Refs. flC|, [lj]], the discrete-symmetry modelQis also AF, in the RG 
language, this means that the added term is an irrelevant perturbation. Let us perform 
a standard RG calculation around the theory with h n = 0. If G^(k\, . . . ,k p ;/3, h n ) is 
the connected p-point correlation function, we can perform an expansion in powers of 
h n , i.e. rewrite 



G (p) (fci, ...,k p ;f3, h n ) = ^KG^ih, . . . , V, 0, . . ■ , 0; A 0), 

g=0 Q - 



(6) 



where G^ p ' q \ki, . . . , kp\ l\, . . . , l q ; (3, 0) is the correlation function with p fields and q 
insertions of the breaking operator at momenta l\, l q , computed for h n = 0, 
i.e. in the 0(3)-model. Such correlation functions cannot be computed directly in 
perturbation theory for l\ = . . . = l q = 0, since the momenta of the insertions are zero 
and perturbative expansions^ are valid only for large momenta, i.e. for distances much 
smaller than the correlation length. Indeed, we are expanding around an ordered 
configuration and this is correct only for |x| <C ^ because of the Mermin- Wagner 
theorem. However, within the usual theoretical framework, there is a standard way out. 
Consider G( p ' q \ki, . . . , k p ; 1%, . . . , l q ; f3, 0) in the perturbative regime (large momenta) 
and use perturbation theory to derive the RG equation 



GM(* 1 ,...,fc„;/ 1 ,...,i„;A0)=0, (7) 



2 In Eq. (|^) we have defined I n so that each polynomial has maximal value one on the vertices of the 
tetrahedron, cube, and icosahedron respectively. The minimal value (maximal in absolute value) of I4 on 
the vertices of the octahedron is —3/2 and of Iq on the vertices of the dodecahedron is —5/9. Of course, the 
normalization of I n is irrelevant for our discussion. 



3 Note that, if a discrete model is obtained fo h n — +00, the dual model (see App. A.l for the definition 
of duality) is obtained for h n = —00. 

4 Strictly speaking, the claim has been made only for \h n \ — 00. However, it is difficult to imagine a 
scenario in which the theory is not 0(3) invariant for h n small and recovers the continuous invariance by 
increasing the strength of the discrete-symmetry term. 

5 Even for large momenta, perturbative expansions are well-defined only for 0(3)-invariant quantities. 
Thus, in order to have an infrared-finite perturbation theory one must rewrite each correlation function in 
terms of 0(3)-invariant expressions. For instance, to compute correlations of we can use the identity 
(Yi 1:mi (cro)Yi 2 , m2 (cr x )) = fyi,j 2 <5mi,m 2 (-Pi:i(fo • (t x ))/(Att), where Pi(x) is a Legendre polynomial. The right- 
hand side is infrared-finite. 
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where t = 1/(3, W(t), -f( n \t), and j(t) are respectively the lattice /3-function and the 
lattice anomalous dimensions of I n and of the field. The crucial assumption is that 
this equation is valid for all momenta, even outside the strictly perturbative region^] 
Taking into account the scaling dimension of the correlation function, at zero ex- 
ternal momenta (ki = Zj = 0), we have 



G^ (0; 0; (3, 0) = A^G^ (0; /?, 0) 



exp 



2- 7 ( n )(s) 



to 



W(s) 



ds 



(8) 



where A^ p,q ^ is a non-perturbative constant. Therefore, we obtain finally 



(0; 0, h n ) = (0; (3, 0) £ -A^f^ exp 



<?=o 



t 2-7< n )(s) 



to 



W{s) 



ds 



(9) 



By using the perturbative RG, we have been able to factor out the h dependence of 
the correlation function. Now, for t — > we find immediately 



GW(0;AM~GW(0;A0)^^^ (h n t^ exp 



9=0 



4ir 



(10) 



0, 



where p n is an easily computable exponent. The correction term diverges for t 
showing that the breaking term is a relevant interaction in the RG sense. 

Equivalently, one can imagine of considering a scale-dependent renormalized cou- 
pling h n {s). Then, the RG flow has the form 



1 dh n (s) 
h n (s) ds 



2-7 {n) [t(s)j 



(11) 



where exp(— s) is the change of the scale and t(s) is the running coupling constant. Since 
-y( n ) (t) vanishes for t — > 0, any perturbation of this type is relevant in the continuum 
limit. We would like to point out that this is not unexpected. Since in two dimensions 
the field a is dimensionless, any polynomial in a is a relevant operator. 

Equations (^) and (|l^) deserve some additional comments. First of all, they give 
an expansion of G^ p '(0; (3, h n ) in powers of the scaling variable 



h n t pn exp(4vr/i) 



(12) 



and are therefore valid only for z <C 1. However, for our purposes, the only relevant 
information we obtain from Eq. ( [To| ) is that the correct scaling variable is z, i.e., that 
we can define the limit t — > 0, h n — ► at fixed z, obtaining 



where <&( p \z) is a nonperturbative crossover function. Eq. (|l 
from the RG equation 

d 



(13) 

can be derived directly 



in) 



(t)h ri 



dh T 



+ |7(*) 



G^\ki, ... ,k n ;f3, h r , 



0, 



(14) 



6 The correlation functions are computable in perturbation theory if we introduce an infrared cutoff, for 
instance if we work in a finite volume. However, they will be infrared divergent, so that one cannot take the 
infinite-volume limit naively. In the standard theoretical framework, one uses the perturbative expressions 
to derive an RG equations and then assumes the such equations are also satisfied by the infinite-volume 
quantities. Whatever the procedure, the result is identical. 



5 



which we assume to be valid for all momenta in a neighborhood of h n = t = 0. Of 
course, we can only compute the scaling behavior of the correlation functions, but not 
their explicit expressions. Our ignorance of the long-wavelength physics is encoded 
here in the nonperturbative nature of the function &( p \ 



Our result (13) provides what in statistical mechanics is called a crossover scaling 
function for a fixed point perturbed by two relevant interactions (or, in this case, by 
one relevant interaction h n and one marginally relevant interaction t): namely, it gives 
the leading behavior in the limit t — > 0, h n — > with z fixed. Note that Eq. (|l~3| ) can 
also be used to predict the behavior of the phase-transition line for h n small ]l5], |i~6| , 
if such a transition exists.[| Indeed, if (3 c {h n ) is the critical point of the theory at fixed 
h n , for h n small we should have, for h n small, 

h n = z*(3 c {h n y n exp[-47r/? c (/ ln )], (15) 

where z* is a nonperturbative constant. 

We want now to explain the numerical results of Refs. fi~0| , [Tl| , ||, who found that 
the difference in behavior between the standard 0(3) model and the discrete model was 
decreasing as (3 increased. A possible explanation of this phenomenon is that the RG 
flux first reduces the size of the perturbation which then increases again as increases. 
Since in the high-temperature regime, one expects indeed the two models to be quite 
similar, this could explain the fact that they are numerically indistinguishable at the 
values of (3 that can be simulated today. 

To make this picture more quantitative, let us consider the Hamiltonian (||) with h n 
and z small, so that we can use Eq. (^). Then, suppose that there exists i eff = l/(3 eS 
such that 7 ( n )(i eff ) = 2 and 2 - y( n \t cS ) < for (3 < f3 eS . In this case, Eq. (|) would 
predict the following behavior. For (3 small, the difference would apparently decrease 
as (3 increases, which could seem to indicate that the interaction is irrelevant. However, 
as soon as (3 becomes larger than f3 the discrepancy starts increasing again. Now, 
notice that the discrete model, in the vicinity of the 0(3) fixed point, will generate 
perturbations of arbitrary spin. However, increases with n, and thus the most rel- 
evant perturbation is associated to I n with the smallest possible value of n. Therefore, 
we should consider n = 3, 4, 6 for the tetrahedron, the cube and the octahedron, the 
dodecahedron and the icosahedron respectively. 

We can try to evaluate /3 eff by using the perturbative expressions for the anoma- 
lous dimension of I n in the 0(3) nearest-neighbor lattice model. Explicit three-loop 
expressions are reported in App. [b]. We obtain the following estimates 



0(3), eff 


= 0.75, 


(16) 


0(4), eff 


= 1.08, 


(17) 


^(6), eff 


= 1.95, 


(18) 



which are quite stable with respect to the loop order. Thus, we expect that for (3 < 
/?(')' eff the breaking to the corresponding subgroup of 0(3) appears as irrelevant. Now, 



7 Eq. (13) is valid for any perturbation and any n, even if it is not associated to a discrete group. For 
instance, one can take n = 1 and Ii((t) = cos 9, obtaining the usual magnetic perturbation. In this case, Eq. 
( |j~3| ) is the usual crossover equation (for p = 1 it is usually called equation of state) reported in textbooks, 
see, e.g., Ref. [|l7|. Instead, Eq. (|l^) that will be derived below makes no sense since at fixed h\ there is no 
phase transition. 
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(3 ~ 1 corresponds to a very small correlation length. Therefore, as soon as £ > 1, one 
immediately realizes that the tetrahedron and the cubic model are not asymptotically 
free. On the other hand, (3 = 1.95 corresponds to a quite large value of the correlation 
length £ = 166.7(4) fl2f . Therefore, we expect the discrepancy to decrease steadily 
as (3 increases, till values of £ of ordeiQ 200 and increase steadily afterwards. Thus, 
a clear signal of the difference between the two models requires simulations on quite 
large lattices with £ S> 200. 

It is important to notice that these estimates are valid for perturbations of the 
lattice nearest-neighbor 0(3) action only. Indeed, as it should be expected, /3^ n )' eff is 
nonuniversal, being the solution of j^(t eS ) = 2, where 7^ n )(t) is the action-dependent 
anomalous dimension of I n . This means that, by changing the action, it is well possible 
that the system shows a non-0(3)-invariant behavior for much smaller values of £. 

The argument given here applies quantitatively only for h n small. Nonetheless, 
for the discrete models considered in the simulations it represents a plausible scenario 
which explains the numerical results and is compatible with the standard theoretical 
framework used in the analyses of the critical behavior of the 0(3) model. Indeed, 
since all simulations are performed in the region in which 1 < £ < 100, it predicts that 
the cubic and the tetrahedron model are clearly different from the 0(3) model, while 
the icosahedron results should mimic the 0(3) ones. 

It is interesting to notice that the values of /3( n )> eff are close to the critical value for 

(n) 

each discrete model 0c ■ For instance, since the tetrahedron model is equivalent to 
the 4-state Potts model with (3 = 3/4/3 Potts we have 

0® = l log 3«0.82, (19) 



which is only slightly higher than (16). Analogously, the cubic model is equivalent to 



the product of three Ising models with (3 = 3/3 Isillg so that 

/3( 4 ) = ^log(l + V2)« 1.32. (20) 

There is also [18] a numerical estimate for (3^ ~ 2.15. Note also that /3( n )' eff < (3^, 
an inequality which shows that our /?( n )' eff correspond always to temperatures above 
the freezing transition. 
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A Discrete subgroups of 0(3) 



8 This number is simply an order of magnitude. Indeed, even if our estimates are quite stable with respect 
to (3, the correlation length varies rapidly with (3. For instance, for /? = 1.90 (resp. 2.00), which are very 
close to our estimate of (3^' cS we have £ = 122.3(3) (228.5(7) resp.) JT§. 
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A.l The Platonic solids 



The icosahedron and the dodecahedron are two Platonic solids. They are regular 



convex polyhedra |TI|, with regular and equal faces, and are such that each vertex 
belongs to the same number of edges. If its faces are p-gons (polygons with p sides), 
q of them surrounding each vertex, the polyhedron is denoted by {p, q}. The possible 
values for p and q may be enumerated as follows. The solid angle at a vertex has q 
face-angles, each (p — 2)n/p. Of course, the sum of these q angles must be less than 
2-7T. Therefore, we have 

111 

- + ->«• 21 
q p 2 

Thus, {p, q} cannot have any other values than 

{3, 3}, {3, 4}, {4, 3}, {3, 5}, {5, 3}. 

They correspond to the tetrahedron, the octahedron, the cube, the icosahedron, and 
the dodecahedron. 

Consider the regular polyhedron {p, q} with its No vertices, Aq edges, and A2 faces, 
where Nq — Aq + A^ = 2 by Euler's formula. As each face touches p edges and each 
edge belongs to 2 faces, then 

pN 2 = 2Aq. (22) 

Analogously, since each vertex belongs to q edges and each edge touches 2 vertices, we 
obtain 

qN = 2Aq. (23) 
These relations and Euler's formula imply 

1111 

- q + - p -2 = N[- < 24 > 



which offers a quantitative evaluation to the inequality fl2T|). 

Consider now the sphere which touches all the edges. If we replace each edge by a 
perpendicular line touching the sphere at the same point, we obtain the Aq edges of 
the dual polyhedron {q,p} which has N2 vertices and A^ faces. 

We are interested in the rotation groups of the regular polyhedra. They are finite 
groups, so that every rotation must have an angle commensurable with ir. In fact, the 
smallest angle of rotations around a given axis is a submultiple of 2tt, and all other 
angles of rotation about the same axis are multiples of the smallest one. If 2n/k is the 
smallest angle, then the rotations about this axis form a cyclic group of order k, and 
one speaks of an axis of k-fold rotation. 

Two dual polyhedra have the same rotation group. The center of the polyhedron 
{p, q} is joined to the vertices, mid-edge points, and centers of faces, by axes of g-fold, 
2-fold, and p-fold rotation. But the vertices, mid-edge points, and centers of faces 
occurs in antipodal pairs. Hence, the total number of rotations, excluding the identity, 
is 

i [(q - 1)N + N 1 + (p- 1)N 2 ] = 2Aq - 1, 

so that the order of the rotation group is 2Aq . 

The rotation group of the icosahedron (and therefore of its dual polyhedron, the 
dodecahedron) has 60 elements. 
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A. 2 The invariants under discrete subgroups 



In this Section we will compute the lowest-degree homogeneous polynomials in the 
(x, y, z) coordinates that are invariant under the action of several discrete subgroups of 
0(3). We will consider only the subgroups which are really three-dimensional isome- 
tries: they are related to the Platonic solids discussed above. The other subgroups of 
proper rotations are the cyclic groups C n , for n > 1, and the dihedric groups D n , for 
n > 2, which are also subgroups of SO (2). 

All invariants of 0(3) can be obtained as powers of the basic degree-two invariant 

h = x 2 + y 2 + z 2 . (25) 



A. 2.1 The cubic group 

Let us fix the cube with vertices at the points (±1, ±1, ±l)/\/3. There are 48 matrices 
which leave invariant the cube, that are 





Twenty- four matrices {2N\ = 24) are proper rotations, while the other 24 matrices 
are obtained by compositions of proper rotations with the antipodal transformation 
diag(— 1, —1,-1). Algebraically, the cubic group is S4 <g) Z2, where S4 is the group of 
permutations of 4 elements. As the cube is dual to the octahedron the two groups of 
invariance are the same. 

It is easy to see that the lowest-order non-trivial polynomial is 

h = x 4 + y 4 + z 4 + all (26) 

On the unit sphere the polynomials can be decomposed into irreducible representations 
of 0(3), i.e. in terms of spherical harmonics. Then 

I *=(l + a ) ^ y o,o + ^v& 4 ,o + \\[^ (*m + n,- 4 ) • (27) 
With the choice a = —3/5, we obtain an operator that is renormalized multiplicatively. 



A. 2. 2 The group of the tetrahedron 

Let us choose the tetrahedron with vertices at the points 

-75(1,1,1), -L(i,-i,-i),-L ( -i,i,-i), ^(-i,-M). 

The matrices which leave invariant the tetrahedron are the 24 matrices of the cubic 
group which have an even number of —1. Remark that there is not the antipodal 
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transformation, but there are 12 proper rotations and 12 reflections. The lowest non- 
trivial polynomial is 



Pu = 


(0,0,1), 


Pd = 


(0,0,-1). 


Pk = 


f 2 

—= cos 
\V5 



I 3 = xyz = -i^ (Y 3 , 2 - y 3 ,_ 2 ) . (28) 
Algebraically this group is A4, the group of even permutations of 4 elements. 

A. 2. 3 The group of the icosahedron 

We shall parametrize the 12 vertices of the icosahedron as follows: 

(29) 

nk\ 2 . fnk\ 1 , 

xJ'7! sm UJ'7! cos( "" 

for k = 1, ■ ■ ■ , 10. In order to construct invariant polynomials under the rotation group 
of the icosahedron, we first consider the cyclic group of order 5 of rotations of 2ir/5 
around the z-axis. On the vertices it acts as a permutation of the form 

(P u ){PiP 3 PbPrP9){PbPAP6PsPio)(Pd) ■ 

In order to determine the invariants under this cyclic group, it is convenient to use 
cylindric coordinates (z, p, 0) € R x R + x [0, 2tt] so that the action of the generator is 

(z, p, 4>) -> (^z, p, <j) + ^ ^ , (30) 

and thus the invariants are z, p and 50. The lowest-order non-trivial polynomial which 
also respect the invariance under the antipodal transformation is of sixth degree and 
has the general form 

I 6 = z (z 5 + az 3 p 2 + bzp 4 + cp 5 cos(50) + dp 5 sin(50)) + elf, (31) 

where we have fixed to one the coefficient of z 6 . Of course, I§ must be the same on all 
the vertices of the icosahedron: this gives the condition 

a = 31 - 46 - 8c. (32) 

We shall then use the cyclic group of order 3 of rotations around the axis which joins 
the origin with the center of the face P u P2Pi- On the vertices it acts as a permutation: 

(.PuP2P4)(PlP 5 P8)(P3P6Pl0)(P 7 P9Pd) ■ 

It gives the conditions 

b = 5, (33) 
c = 2, (34) 
d = 0, (35) 

so that a = — 5. We obtain finally 

h = z (z 5 - 5z 3 p 2 + 5zp 4 + 2p 5 cos(50)) + e (z 2 + p 2 f (36) 
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or, in Cartesian coordinates, 



(z 5 - 5z 3 (x 2 + y 2 ) + 5z(x 2 + y 2 ) 2 + 2x(x 4 - Wx 2 y 2 + by 4 )) 



+e (x 2 + y 2 + z 2 ) 3 



(37) 



It can be checked that the all the other transformations, like the cyclic groups of 
order 2 around the center of an edge, for example the transformation which acts as a 
permutation like 

(P U P 2 ) (PxP 6 ) (P 3 P 8 ) (P 4 P W ) (P 5 P 9 ) (P 7 P d ), 

leave Iq invariant. Also the antipodal transformation leaves Iq invariant. 
On the unit sphere, in terms of spherical coordinates, we obtain 



21 



16 Mvr 



h=\— + e\ V47ry ,o + — \l — *6.o - — 



16 4vr 



21 V 13 



1001 



(y 6)5 - y 6 _ 5 ) 



Therefore, for 



5 

21 



(38) 



(39) 



we obtain an operator which is multiplicatively renormalized. 

Of course we could use a different position of the icosahedron in space, for example 
we could take for the vertices 



1 



V2T 



(0,±r,±l) 



1 



(±T,±1,0) 



1 



(±1,0, ±r) 



(40) 



where r is golden ratio. In this basis the cyclic group of rotations of order five is 
generated by 

r t-1 -1 \ 
r - 1 1 r . (41) 
1 -r t-1 



The new choice can be recovered from the old one by a rotation and this produces a 
different polynomial which turns out to be 



I'* 



21 / ' V 273 



1 



1 



6,0 



~ (y 6 , 4 + y 6 ,-4) + ^ V IT (y6 ' 6 + y6 -- 6) 



(42) 



which, of course, with the choice (^9|) belongs to the same multiplet I = 6. 

The icosahedron is dual to the dodecahedron and thus their group of invariance is 
the same. 



B The perturbative results 

In j2(| the anomalous dimension of all non-derivative dimension-zero operators was 
computed for the nearest-neighbor lattice 0(N) cr-model up to three loops. For generic 
N, a suitable basis for these operators is given by 

, = a jl . . . <r jn - traces, (43) 



11 



(n) 

where "traces" must be such that - n is completely symmetric and traceless. These 
polynomials are irreducible 0(./V)-tensors of rank n and as such they renormalize mul- 
tiplicatively with no off-diagonal mixing. For N = 3 this representation is equivalent 
to that of the spherical harmonics. 

The anomalous dimension 7 ( n ) (i) of these operators is given by 



7 



(n) 



(t) = t ] t + it ] t 2 + ii n) t 3 + 0(i 4 ), 



(44) 



where 



where 



(n) 

7o 

(n) 

7i 

(n) 
72 



n(iV + n-2) 

4^ ' 

n(iV + n- 2) 

167T ^ 

n(N + n-2) 

47T 



I r d d p d d q 
1_ ~4 ./ (2vr) d (2vr) rf 



iV 



4vr 2 



^ 4 - T 



(P + <?) 5 



+ 



11 

96 



~2 ~2 

p — q 



p 2 q 2 [(p + (?) 2 ] 2 



(45) 



(46) 



and 



Numerically G\ 



p 2 = EpI = £( 2sin f) 2 - 

: 0.0461636 . For N = 3 we obtain 

7 ( 3 )(t) = 0.95493t + 0.238732t 2 +0.135755t 3 , 
7 ( 4 )(i) = 1.59155* + 0.397887t 2 + 0.226258* 3 , 
7 ( 6 ) (t) = 3.34225t + 0.835563t 2 + 0.475142i 3 . 



(47) 



(48) 
(49) 
(50) 
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